Chiral exponents in frustrated spin models with noncollinear order 
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We compute the chiral critical exponents for the chiral transition in frustrated two- and three- 
component spin systems with noncollinear order, such as stacked triangular antiferromagnets 
(STA). For this purpose, we calculate and analyze the six-loop field-theoretical expansion of the 
renormalization-group function associated with the chiral operator. The results are in satisfactory 
agreement with those obtained in the recent experiment on the XY STA CsMnBra reported by V. 
P. Plakhty et al., Phys. Rev. Lett. 85, 3942 (2000), providing further support for the continuous 
nature of the chiral transition. 
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The critical behavior of frustrated spin systems with 
noncollinear order is still a controversial issue, field- 
theoretical (FT) methods, Monte Carlo (MC) simula- 
tions, and experiments providing contradictory results 
in many cases. At present there is no agreement on the 
nature of the phase transition, and in particular on the 
existence of a new chiral universality class See, e.g., 
the recent works and Refs. for reviews. 

In magnets noncollinear order is due to frustration that 
may arise either because of the special geometry of the 
lattice, or from the competition of different kinds of inter- 
actions. Typical examples of systems of the first type are 
two- and three-component antiferromagnets on stacked 
triangular lattices Their behavior at the chiral tran- 
sition may be modeled by a short-ranged Hamiltonian 
for iV-component spin variables Sa, defined on a stacked 
triangular lattice as 

HSTA = ^(*) ■ ^(^■) " E ^(*) ■ ^(^■)' 

where J < 0, the first sum is over nearest-neighbor pairs 
within triangular layers, and the second one is over or- 
thogonal interlayer nearest neighbors. Frustration due to 
the competition of interactions is realized in helimagnets. 

In these models frustration is partially released by mu- 
tual spin canting and the degeneracy of the ground state 
is limited to global 0(iV) spin rotations and reflections. 
At criticality one expects a breakdown of the symmetry 
from 0{N) in the high-temperature phase to 0{N — 2) in 
the low-temperature phase, implying a matrix-like order 
parameter. In particular, the ground state of the XY 
systems shows the 120° structure of Fig. |l|, and it is Z2 
chirally degenerate according to whether the noncollinear 
spin configuration is right- or left-handed. The chiral de- 
grees of freedom are related to the local quantity [|] 

Cab OC ^ Sa(t)Sb{3) ~ Sb{i)Sa{3) (2) 
A 

where the summation runs over the three bonds of the 
given triangle. The definition of Cab can be straightfor- 
wardly generalized to the case of A'^-component spins. 




FIG. 1. The ground-state configuration of three XY spins 
on a triangle coupled antiferromagnetically. 

Many experiments are consistent with a second-order 
phase transition belonging to a new (chiral) universality 
class (see, e.g., Refs. [^,|| for reviews). Further experi- 
mental evidence in favor of a chiral continuous transition 
has been recently reported in Ref. |^ , showing the simul- 
taneous occurrence of spin and chiral order in the XY 
stacked triangular antiferromagnet (STA) CsMnBrs. 

On the theoretical side the issue has been controver- 
sial. MC simulations ||9|-p^ (see Refs. [^,|^ for reviews) 
have not been conclusive in setting the question. Simula- 
tions of STA's are consistent with continuous transitions, 
but with critical exponents that are not in a satisfactory 
quantitative agreement. In Ref. jl^ the results for the 
XY STA are interpreted as an evidence for mean-field 
tricritical behavior. Moreover, MC investigations [|l^ of 
special lattice spin systems, that on the basis of their 
symmetry should belong to the chiral universality class, 
show clearly a first-order transition. 

In a recent Letter 

i the issue has been studied by 
a continuous renormalization-group (RG) approach (see 
also Refs. [^|7|). The results favor a first-order transi- 
tion, since no evidence of stable fixed points is found. Ac- 
cording to this first-order transition picture, the apparent 
continuous critical phenomena observed in experiments 
are interpreted as first-order transitions, weak enough 
to effectively appear as second-order ones. Note however 
that the practical implementation of this method requires 
an approximation and/or truncations of the effective ac- 
tion. So these studies may not be conclusive. 

FT studies of systems with noncollinear order are 
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FIG. 2. RG flow in the quartic couplings {u, v) plane for 
A'^ = 2, 3. It shows the stable chiral fixed point denoted by C, 
and the unstable antichiral {A), 0{2N) Heisenberg {H) and 
Gaussian (G) ones. 

based on the Landau-Ginzburg- Wilson 0(A'^)xO(2)- 
symmetric Hamiltonian [QH] 

a a 
a, 6 

where (pa (1 < a < 2) are two sets of TV-component vec- 
tors. Frustrated XY and Heisenberg spin systems with 
noncoUinear ordering, such as STA's, are described re- 
spectively by the N = 2 and = 3 case with vq > 0. 
The presence of a stable chiral fixed point, conjectured by 
Kawamura jl],^ , has been recently confirmed by the anal- 
ysis of the perturbative six-loop series in the framework 
of the fixed-dimension expansion . As sketched in 
Fig. ^, a stable chiral fixed point C appears for both 
XY and Heisenberg cases. The critical exponents char- 
acterizing the stable chiral fixed point turn out to be in 
satisfactory agreement with experiments. Note that in 
this RG picture first-order transitions are still possible 
for systems that are outside the attraction domain of the 
chiral fixed point. In this case, the RG flow runs away 
to a first-order transition. This may explain some exper- 
iments (for example those for the CsCuCls compound, 
see, e.g., Refs. 1^,^) and MC studies for special lattice 
systems [p^, where first-order transitions are observed. 

Beside the conventional critical exponents /3, 7, i^, etc... 
related to the standard spin order, one may consider ad- 
ditional critical exponents related to the behavior of the 
chiral degrees of freedom. If spin and chiral order oc- 
cur simultaneously, one expects 1/^ = 1^ where Vc is the 
exponent associated with the correlation length defined 
from the chiral correlation function. Introducing a chiral 



external field he coupled with the chirality Cabi one may 
write the singular part of the free energy as |l) 

F,ins^t^-"f{h/t^,hjt^^), (4) 

where t is the reduced temperature, A — f3 + j, and pc is 
the chiral crossover exponent. Then, differentiating with 
respect to one may obtain the RG relations 

f3c=^3iy - (j)c, 7c = 20c - 3j/, (5) 

where Pc and 7c describe respectively the critical behav- 
ior of the average chirality and of the chiral susceptibility. 

The first esperimental estimates of the chiral exponents 
(pc and have been recently reported in Ref. |^ for the 
transition of the XY STA CsMnBrg: 

0c - 1.28(7), /3c = 0.44(2), (6) 

measured respectively in the high- and low-temperature 
phase. On the theoretical side, there are a few MC results 
for the STA spin models (P, and very little from field- 
theoretical approaches. The chiral exponents have been 
only computed to 0(1/7V) and 0(e) in the correspond- 
ing expansion frameworks However, these results do 
not allow a quantitative comparison, essentially for two 
reasons: because the series are too short and, most im- 
portantly, as discussed in Ref. |^ , the chiral fixed point 
for the XY and Heisenberg cases is not analytically con- 
nected with the one found in the large-iV and small-e 
region. In order to obtain results that can be compared 
with experiments, one should compute them directly for 
d = 3 and for the number of components of interest, i.e. 
TV = 2,3. 

In this Letter we compute the chiral exponents using 
the FT approach of Ref. Q|, i.e. by computing and an- 
alyzing the six-loop perturbative expansion of the chiral 
RG functions. In the fixed-dimension FT approach one 
performs an expansion in powers of appropriately defined 
zero-momentum quartic couplings (see, e.g., Ref. § and 
references therein). In order to obtain estimates of the 
universal critical quantities, the perturbative series are 
resummed and then evaluated at the fixed-point values 
of the couplings. The comparison with the experimental 
results (|^) will represent a highly nontrivial check of the 
FT description of the transition and of the Kawamura's 
conjecture that these systems undergo continuous tran- 
sitions belonging to distinct chiral universality classes. 

In order to compute the universal quantities character- 
izing the critical behavior in the high-temperature phase, 
one introduces a set of zero-momentum conditions for the 
(one-particle irreducible) two-point and four-point corre- 
lation functions (see, e.g., Ref. for details), which re- 
late the zero-momentum quartic couplings u and v and 
the mass scale m to the corresponding Hamiltonian pa- 
rameters mqj ^^0 and r. In particular, 

(P) = 5a,6,jZ^' W + + 0(/)] . (7) 
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In addition, one defines the function Zt througli the re- 
lation r^^^;^](0) = bahkjZ^^, where T^i'^) jg the (one- 
particle irreducible) two-point function with an insertion 
of \<iP' ■ The fixed points of the theory are given by the 
common zeros m*, v* of the /3-functions 

Pu(u,v) = m——, iy(u,v)=m——, (8) 
dm am 

calculated keeping uq and vq fixed. The critical behavior 
is determined by the stable fixed point of the theory. 
The analysis of the six-loop expansion of the /3-functions 
provided a rather robust evidence of the existence of a 
stable fixed point as shown in Fig. ||. The critical 
exponents 77 and i' are then derived by evaluating the 
RG functions 



dlnm 



rit{u,v) 



dlnZt 
dlnm 



(9) 



at the chiral fixed point u*, v*. The resulting exponents 
are v = 0.57(3), r] = 0.09(1), 7 = 1.13(5) for the XY 
case, and v = 0.55(3), 77 = 0.10(1), 7 = 1.06(5) for the 
Heisenberg case j2^] , which are in substantial agreement 
with the experimental results. 

In order to evaluate the chiral exponents, we consider 
the operator 

Cckdl{x) = 4'ck{x)(t)di{x) - (f>ci{x)(f>dk{x), (10) 

and define a related renormalization function Zc from the 
one-particle irreducible two-point function F^^'^^ with an 
insertion of the operator Cai,bj, i-e. 

r^^'^^(0)aj,bj,cfe(ii = Z^^ Tabcd,ijkl 

where 



(11) 



Tabcd,ijkl — {SacSbd ^ Sadhc) {^ik^jl — Sil&jk) (12) 

SO that Zc(0,0) — 1. Then, we compute the RG function 
d\nZr 



dim 



du 



dv 



and its value rjc s^t u — u* , v — v* , where u*, v* is the 
position of the stable chiral fixed point Finally, the 
RG scaling relation 



4>c ^ {2 + rjc - ri) v 



(14) 



allows us to determine (pc- 
We computed F^'^'^) (g) 

to six loops. The calculation is 
rather cumbersome, since it requires the evaluation of 563 
Feynman diagrams. We handled it with a symbolic ma- 
nipulation program, which generates the diagrams and 
computes the symmetry and group factors of each of 
them. We used the numerical results compiled in Ref. pl| ] 
for the integrals associated with each diagram. The re- 
summation of the series was performed using the method 



outlined in Refs. [^2|jj]. The very lengthy expression of 
the six-loop expansion of r]c{u,v), details of its calcula- 
tion, and its analysis will be reported elsewhere. The 
results of our analysis are 



1.43(4) for XY, 
1.27(4) for Heisenberg. 



(15) 
(16) 



The errors are indicative of the spread of the results 
yielded by the analysis when changing the resummation 
parameters and varying the location of the chiral fixed 
point within the range reported in Ref. Q . Using the RG 
relations (||) and the estimates of , one may also de- 
rive corresponding results for the other chiral exponents, 
obtaining for example 



Pc = 0.28(10) 
Pc = 0.38(10) 



for XY, 

for Heisenberg. 



(17) 
(18) 



We may compare these results with the experimental 
ones (H) . Our estimate of 4>c is somewhat higher than the 
estimate (^) while the estimate of Pc is correspondingly 
somewhat lower. In any case, the difference is of the or- 
der of one combined error bar. We may also compare our 
results to the available MC estimates for the XY STA 
spin model, that are Pc = 0.45(2), 7c = 0.77(5), 0c = 
1.22(6) from Ref. §, and Pc = 0.38(2), 7c = 0.90(9), 
(j)c — 1.28(10) from Ref. ||ll[ (in this work a mean-field tri- 
critical behavior is conjectured for the transition). These 
results are close to the experimental ones and thus show 
the same deviations with respect to our FT results. 

For TV = 3 we can compare our results with the MC 
ones for the three-component STA spin model, that are 
Pc = 0.55(4), 7c = 0.72(8), (jjc = 1.27(9) from Ref. § 
and Pc = 0.50(2) 7c 0.82(4) and 0c = 1.32(5) from 
Ref. The FT estimate of (pc is in perfect agreement 
with the MC results, while the estimate of Pc is some- 
what lower, although compatible within error bars. Ap- 
parently, this is due to the fact that our estimate of v is 
somewhat lower than those obtained in MC simulations. 

In conclusion, the FT results are in satisfactory agree- 
ment with the experimental and MC estimates. This is a 
nontrivial check of the FT approach, shows its predictive 
power in spite of the fact that the perturbative series are 
not Borel summable — still we take into account the lead- 
ing diverging behavior, see Ref. [Q — and strengthens the 
evidence for the continuous nature of the chiral transition 
in XY and Heisenberg STA's. 
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